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One way to increase the aircraft fuel efficiency is to reduce structural weight while 
maintaining adequate structural airworthiness, both statically and aeroelastically. A design 
process which incorporates the object-oriented multidisciplinary design, analysis, and 
optimization (MDAO) tool and the aeroelastic effects of high fidelity finite element models to 
characterize the design space was successfully developed and established. This paper 
presents two multidisciplinary design optimization studies using an object-oriented MDAO 
tool developed at NASA Armstrong Flight Research Center. The first study demonstrates 
the use of aeroelastic tailoring concepts to minimize the structural weight while meeting the 
design requirements including strength, buckling, and flutter. Such an approach exploits the 
anisotropic capabilities of the fiber composite materials chosen for this analytical exercise 
with ply stacking sequence. A hybrid and discretization optimization approach improves 
accuracy and computational efficiency of a global optimization algorithm. The second study 
presents a flutter mass balancing optimization study for the fabricated flexible wing of the 
X-56A model since a desired flutter speed band is required for the active flutter suppression 
demonstration during flight testing. The results of the second study provide guidance to 
modify the wing design and move the design flutter speeds back into the flight envelope so 
that the original objective of X-56A flight test can be accomplished successfully. The second 
case also demonstrates that the object-oriented MDAO tool can handle multiple analytical 
configurations in a single optimization run. 


Nomenclature 


AFRC = Armstrong Flight Research Center 

AFRL = Air Force Research Laboratory 

ARMD = Aeronautics Research Mission Directorate 
AW IB = anti-symmetric wing first bending 

AW IT = anti-symmetric wing first torsion 

AWBT = anti-symmetric wing bending torsion 

BLF = buckling load factor 

CDV = continuous design variable 

CTMIN = constraint tolerance for inequality constraints in DOT 

CG = center of gravity 

DDV = discrete design variable 

DOT = design optimization tools 

EFEW = Empty Fuel Empty Water 

FFFW = Full Fuel Full Water 

FI = failure index 

GA = genetic algorithm 

GVT = ground vibration test 

KEAS = knots equivalent air speed 

LMSW = Lockheed Martin Skunk Works 

MDAO = multidisciplinary design, analysis, and optimization 
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MSC = MacNeal-Schwendler Corporation 

NASA = National Aeronautics and Space Administration 

SW1B = symmetric wing first bending 

SW1T = symmetric wing first torsion 

SBFF = symmetric body freedom flutter 

SWBT = symmetric wing bending torsion 

V = aircraft velocity 

V-f = velocity-frequency 

V-g = velocity-damping 

V L = aircraft limit speed 

I. Introduction 

T HE modern aircraft design of today is a multidisciplinary design procedure that requires a higher level of 
process integration, design innovations, and revolutionary design concepts. One of the major goals of the Fixed 
Wing and High Speed Projects under the National Aeronautics and Space Administration (NASA) Aeronautics 
Research Mission Directorate (ARMD) is the development of improved prediction methods and technologies for 
lower noise, lower emissions, lower sonic boom, and higher performance for subsonic, supersonic, and hypersonic 
aircraft. One way to increase the fuel efficiency is to reduce structural weight while maintaining adequate structural 
airworthiness, both statically and aeroelastically. To meet these goals, a systematic multidisciplinary design tool will 
support the design of a next generation aircraft. Most commercially available multidisciplinary design, analysis, and 
optimization (MDAO) tools, such as in MSC Nastran (MSC Software Corporation, Newport Beach, California), 
have been developed to perform within a limited number of disciplines with a single fidelity modeling capability and 
a single analytical model. Although a multitude of tools have been developed and well adapted for various 
disciplines in aircraft design/analysis, they have not been fully integrated to work together naturally. Data transfers 
between the codes are performed by using the multidisciplinary computing environment. Several works in 
aeroelastic optimization have been published using academic codes and/or commercial codes with high-fidelity 
and/or tight coupling. 1-3 

Over the past several years, the NASA Armstrong Flight Research Center (AFRC) (Edwards, California) has 
developed an object-oriented MDAO tool. 4 This tool provides a computational environment in which the optimizer 
can effectively receive objective and constraint function values from various disciplines through the interface 
variables. This tool has been used to evaluate several real-world optimization problems including structural dynamic 
model tuning, 5 ’ 6 unsteady aerodynamic model tuning, 7 and mishap investigations. While the MDAO tool is a 
promising technology for complete vehicle system level design and optimization, 8 many gaps and challenges remain 
to be addressed. One of the main challenges is the high computational cost of analysis which is due to the 
demanding nature of the high fidelity analysis in terms of model fidelity and computational resources required, as 
well as organizational complexity. 4 As the level of fidelity is increased, the number of design variables and 
constraints increases. Additionally, the interdisciplinary coupling inherent in MDAO tends to present additional 
challenges beyond those encountered in a single-discipline optimization. 

With the advanced composites materials of today, and developed optimization techniques combined, the MDAO 
tool has become a powerful technology allowing structures to be designed more efficiently to meet structural or 
aircraft performance requirements. Composite optimization with aeroelastic considerations has been a research topic 
for many years. 10 " 13 Composite materials are employed to lower the structural weight due to their high strength-to- 
weight ratio over that of metals. An overview of the earlier research provided by Vanderplaats and Weisshaar 14 and 
Ringertz 15 presents a mass minimization with aeroelastic constraints. In the past, the aeroelastic tailoring theory 12 
combined with the advanced composites in aircraft construction have been successfully applied to improve the 
performance of modern aircraft (e.g. the NASA X-29 aircraft). 16 

The primary objective of the aeroelastic optimization studies in this paper is to demonstrate the capability of the 
MDAO tool for leveraging selected engineering disciplinary codes along with the high fidelity aeroelastic analysis. 
A high fidelity finite element model of the X-56A aircraft was chosen for the optimization studies. This paper 
describes the results of two optimization studies. The first optimization study concentrates on designing a wing to 
avoid flutter within the flight envelope using aeroelastic tailoring concepts to minimize the structural weight while 
meeting the design requirements including strength, buckling, and flutter. This study is demonstrated by using a 
hybrid optimization technique which first uses the global optimizer to find the global optimum values and then 
switches to a gradient-based optimizer to accelerate the optimization procedure. The second optimization study 
restrains the vehicle flutter speeds within a desired range so that during flight testing engineers can examine an 
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active flutter suppression system within the flight envelope. This second study also demonstrated the capability of 
the MDAO tool to include multiple analytical models (multi-point design) in a single optimization run. Multi-point, 
multi-mission optimization is considered a necessity in the aircraft design. 17 This optimization study was performed 
using a gradient-based optimizer and was followed by sensitivity studies. Two different weight configurations were 
used in this study. 

II. Object-Oriented Multidisciplinary Design, Analysis and Optimization Tool 

Over the past several years, the NASA AFRC has developed an object-oriented MDAO tool. This tool is 
designed to handle the complex optimization problems that involve multi-point, multi-mission, multi-level analyses 
and multi-fidelity modeling for multi-disciplinary analyses. 8 As shown in Fig. 1, each of the discipline modules is a 
generic system that can be easily developed, maintained, and extended. This MDAO tool leverages existing in-house 
codes and third party analysis tools such as Nastran, ZAERO (Zona Technology Inc., Scottsdale, Arizona), 18 and a 
CFD solver. Users can choose an optimization methodology and provide initial design values as well as side 
constraints for both continuous and discrete design variables. The MDAO tool allows users easy information access, 
problem formulation, and execution to enable multidisciplinary optimization during the preliminary and detailed 
design stages of subsonic, transonic, and supersonic aircraft. Three optimization algorithms have been incorporated 
into the MDAO tool to include a gradient -based algorithm, a genetic algorithm (GA), 19 and the big-bang big-crunch 
algorithm. 20 In the Object-Oriented MDAO Tool, design optimization tools (DOT) 21 are used for the gradient-based 
algorithm. 

The structural design of an airframe in MDAO is determined by multidisciplinary criteria which include 
strength, stiffness, and buckling loads for sizing optimization under given loading conditions and design 
requirements. Structural analysis is a complex effort involving trim analysis, loads analysis, modal analysis, and 
aeroelastic as well as aeroservoelastic analyses in generating output such as structural weight, deformations, 
divergence speed, open and closed loop flutter speeds and frequencies, natural frequencies, mode shapes, and 
gain/phase margins. 

Several sub-modules for structural analysis disciplines have been developed and implemented in the MDAO 
environment. These sub-modules automate data calculation, provide data output in useful formats, and provide 
information to compute the objective and constraint function values. Although these sub-modules are developed 
mainly for MSC Nastran 22 and ZAERO, they can be easily customized for other analytical tools. Some of these 
modules are discussed below. 

1) Weight Sub-module-MSC Nastran is used to determine the model total weight, moment of inertias, and 
location of the center-of-gravity (CG). 

2) Strength Sub-module-The material stress constraint imposed in the design process is based upon the 
material strength design allowable. For a feasible design, all stresses in the structure must be less than the 
design allowable stress of the material for all load cases. MSC Nastran was used to compute the element 
stress, strain, tip deflection, and composite ply failure indices. Failure theories are calculated at the ply 
level and the result for each element and for each ply is retrieved from the output file. 

3) Buckling Sub-module-MSC Nastran is used to calculate the eigenvalue associated with the first buckling 
mode of the structural design with the given static load. The critical buckling load factor can also be 
calculated for the structure. 

4) Modal Analysis Sub-module-MSC Nastran is used to determine the modal characteristics (natural 
frequencies and mode shapes) of the structural designs. These data are then passed on for flutter analysis 
with frequency constraints if needed. 

5) Divergence and Flutter Sub-module-These analyses are performed using either ZONA’s ZAERO code 
with g-method 18 solution or MSC Nastran code with PK solution method 23 for computing the divergence 
and flutter speeds. These speeds are retrieved from the output file for each configuration and flight 
condition. The final flutter results are displayed as velocity-damping (V-g) and velocity-frequency (V-f) 
plots and used to examine the flutter mechanism. The flow chart of the flutter analysis module in the 
MDAO tool is presented in Fig. 1. 

III. Aeroelastic Model 

The X-56A aircraft is an unmanned aircraft being developed by the U.S. Air Force Research Laboratory (AFRL) 
(Wright-Patterson Air Force Base, Ohio) to test technologies that will be needed for new kinds of lightweight, 
flexible aircraft. The 7.5 foot-long aircraft has a 28-foot wing span and is powered by two turbine engines. The 
wings have an aspect ratio of 14. The aircraft is being built and will be flown by Lockheed Martin Skunk Works® 
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(LMSW) (Lockheed Martin Corporation, Bethesda, Maryland) under contract with the AFRL. After these first flight 
tests are completed, the valuable AFRL assets, two reusable center bodies, one stiff wing, three flex wings, and 
ground control stations as shown in Fig. 2, will stay at NASA AFRC for future technology demonstrations. The 
initial purpose of the X-56A aircraft is to demonstrate the simultaneous active suppression of three flutter 
mechanisms: first symmetric body- freedom flutter (SBFF), first symmetric wing bending-torsion (SWBT) flutter, 
and first anti-symmetric wing bending-torsion (AWBT) flutter, through the use of feedback controls. Thus, this 
aircraft is designed to have three flutter speeds within its flight envelope. Further details of the aircraft design can be 
found in reference. 24 

The X-56A structural model was generated using the MSC Nastran code. The structural finite element model is 
shown in Fig. 3. There are two straight spars and 16 straight ribs for each wing half. The wing skin, spars, and ribs 
are modeled using shell elements with composite material properties. The fuel and ballast weights are modeled using 
concentrated mass elements connected to the ribs and spars by multiple point constraint elements. The wing is 
connected to the center body by point spring elements. The ply orientation is defined in the local material orientation 
coordinate system. Positive angle direction is defined using the right-hand rule about the local z-axis. 

In this study, the unsteady aerodynamics were modeled using the Doublet Lattice method IN in MSC Nastran. 
This model produces a matrix of linear aerodynamic influence coefficients that describe the unsteady pressure 
changes of the aerodynamic degrees-of-freedom. Flutter analyses is performed using the MSC Nastran PK solution 
method which interconnects the structural and aerodynamic grids by spline interpolation and finds the generalized 
aerodynamic matrix using the structural modal matrix. 

IV. Description of the Optimization Problem Formulation 

The constrained optimization problem is defined as follows in Eqs. (1)— (4): 


Find the design variables X = 



that minimize the objective function, 

/(X) = total structural weight or flutter speeds (1) 


such that: 

inequality constraints, 


equality constraints, 


and side constraints, 


5j (X) <0, j=l,M 
h k (X) = 0, k = l,L 
x\ <x t < xj 7 , i = 1,1V 


where xf and xf are the lower and upper bounds on each of the design variables, respectively. 


( 2 ) 

( 3 ) 

( 4 ) 


A. Design Variables 

Design variables are used in the iterative design process to make changes in structural sizes, properties, and 
geometry. The number and type of design variables depend upon the outer shape of the aircraft as well as the 
geometry parameterization of the spars, ribs, upper and lower wing skin, and internal wing structure topology. 


B. Objective, Response, and Constraint Functions 

The response functions considered can be divided into structural and aeroelastic responses. Each of these 
responses can be considered as an objective or a constraint for the problem. In general, the structural responses 
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include structural weight, von Mises stress and strain distributions, displacement, buckling loads factor, modal 
frequencies, mode shapes, et cetera. The aeroelastic responses include divergence speed, flutter speeds, and flutter 
frequencies. The number and type of response functions to be considered depend upon the design requirements. The 
behavior constraints could be equality constraints or inequality constraints, depending on the nature of the problem. 

V. Design Process and Results 

For the first optimization study, the initial baseline model originated with the LMSW X-56A aircraft with flex 
wings finite element model with strengthened wing upper and lower skin thicknesses. The initial model for the 
second optimization study was based on the GVT correlated finite element model (validated model). The structural 
analyses were computed by MSC Nastran codes and the response functions or performance indices were extracted 
from the Nastran output files. The flutter onset was computed by MSC Nastran codes, and the V-g and V-f plots 
were used to examine the flutter mechanism. The V-g and V-f plots indicate visually which natural vibration modes 
are responsible for the aeroelastic coupling. For this study, the point where the aerodynamic damping value crosses 
the 3% line is defined as the flutter onset point to account for an assumed structural damping. The first twenty 
normal modes were used to represent the structure for the flutter analysis. The flutter speed was determined at a 
Mach number of 0.16. 

A. Study 1: Aeroelastic Tailoring 

This study concentrates on designing a wing to avoid flutter onset within the flight envelope (flutter speed 
greater than 1.15 VL) using aeroelastic tailoring concepts to minimize the structural weight while meeting the design 
requirements. The objective function is structural weight. Inequality constraints were used in this study, which 
include composite failure index (FI), buckling loads factor (BLF), flutter speeds, and flutter frequencies. The 
constraints and constraint tolerance values are summarized in Table 1. 

A constraint is violated if its numerical value is more positive than the tolerance value. The flutter constraint 
tolerance for normalized inequality constraint violation was set to 0.002. For DOT optimization, an additional 
constraint tolerance parameter value, CTMIN, 14 of 0.003 was applied. The wing is designed to operate at 2.5 g 
and -1.0 g and a factor of safety of 1.5 maneuver load conditions. These static maneuver loads are used as the 
applied loads for the buckling load factor calculation as well. The flutter speed requirement is that the lowest flutter 
speed must be greater than normalized speed of 1.62 to ensure it is flutter free within the flight envelope. 

For the studies herein, a high fidelity structural model was created based on the LMSW X-56A finite element 
model as described in Section III. The Empty Fuel Empty Water ballast (EFEW) weight configuration was used as 
an initial baseline model. The outer shape of the aircraft was fixed (i.e. the aerodynamic shape did not change), and 
the only design variables were the ones related to the wing skin lamination parameters including ply thickness and 
ply orientation angles. To avoid an unrealistic thickness distribution and reduce the number of design variables in 
the optimization process, where each element of the finite element model is a design variable, the upper and lower 
skin panels were each divided into two design zones, as shown in Fig. 4. There are six design variables including 
three ply thickness and three ply orientation angles in each zone. The thickness and angle of plies are assumed to be 
constant throughout each zone. Also, design variable linking was used for the left and right wings to reduce the 
number of design variables and ensure symmetric design and mass distribution. 

Two optimization cases, 1 and 2, with different design variable sets listed in Table 2 were investigated. The 
optimizations were performed with the minimization of the weight of wing composite laminates consisting of a total 
of 12 design variables (ply thickness only) and 24 design variables (12 ply thickness and 12 ply orientation angles) 
in cases 1 and 2, respectively. A hybrid optimization approach, first using a genetic algorithm (GA) to have a rough 
estimate of the global optimum solution and then using the DOT optimizer to accelerate the computational speed, 
was used in this study. In a GA, the population size depends on desired resolution and the number of design 
variables. In general, several hundreds or thousands of possible solutions are required to cover the entire range of 
possible solutions, which is very time-consuming. Therefore, first the GA optimizer with discrete design variables 
(DDV) was used to locate the optimum region. The GA with a population size of 200 and 30 generations with 
identical convergence criteria were used for cases 1 and 2. Then DOT with continuous design variables (CDV) was 
employed to accelerate the optimization and fine tune the optimum design. The starting design for step 2, DOT 
optimization, is based on the final design from step 1, GA optimization. 

The structural response values and normalized flutter response values of the initial design are listed in Table 3 
and Table 4 for cases 1 and 2, respectively. The normalized total weight of the initial design is set to 1.0. The 
starting design variable values and side constraints are listed in Table 5 and Table 6. The skin thickness upper limit 
was used for all the thickness design variables to guarantee that the optimization starts from a feasible region. Cases 
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1 and 2 have the same flutter mechanism, but differ slightly in flutter speed and frequency. The difference is due to 
the different initial ply angles. The V-g and V-f plots of the initial design are shown in Fig. 5. The SBFF, coupling 
of the short period mode and the symmetric wing first bending mode, were at approximately 1.69 and 1.92 for cases 
1 and 2, respectively. The V-g and V-f plots also show that the coupling of the symmetric wing first bending mode 
and the symmetric wing first torsion mode has led to the SWBT normalized flutter speeds at approximately 2.18. 
Note that, AWBT normalized flutter speeds occurred at approximately 2.38. 

The optimization results of steps 1 and 2 are summarized in Table 3 and Table 4 for cases 1 and 2, respectively. 
As seen in the tables all of the structural response values satisfied the constraints within tolerance. The normalized 
final weights of cases 1 and 2, after global optimization (step 1), are 0.88 and 0.78, respectively. The step 1 results 
show that case 2, which included ply orientation angles as design variables, produce a better design weight than case 
1 (i.e. additional 10% weight reduction). The additional weight saving is mainly due to the aeroelastic tailoring 
effect. In step 2 (DOT optimization), the optimizer is able to lower both cases of the normalized design weight 
further to 0.87 and 0.72. The number of function calls for cases 1 and 2 are 1406 and 3416 in step 1, and 77 and 297 
in step 2, respectively. Overall Step 2 converged much faster than the GA optimization in step 1, and it only 
consumed 5% to 8% of the total computational cost. The final design variable values after step 1 and step 2 for cases 
1 and 2 are listed in Table 5 and Table 6. 

In the airframe industry, the thickness of each type of composite ply is predefined. The final design result from 
DOT optimization with CD Vs is not practical for manufacturing. The ply thickness has to be discretized by either 
rounding up or rounding down to the predefined value, and the fabric ply is placed at a 15-degree discretized angle 
for fabrication. Therefore, the discretization was performed after the DOT optimization. The round up and round 
down design variable values for cases 1 and 2 are listed in Table 5 and Table 6. The response values based on the 
discretization approach are listed in Table 3 and Table 4. The results show that the round down approach for both 
cases leads to a violated design. The critical flutter constraint, which is the active constraint, became violated (i.e. 
1.60 and 1.42 were less than the critical flutter speed requirement of about 1.62). 

For the round up approach, although the critical flutter speed is satisfied for both cases with a weight penalty, 
other response values changed, which could be good or bad (i.e. decrease in buckling load factor response value). 
There is a possibility that the final design will not satisfy all the design requirements after the discretization. For 
case 2, Table 4 shows that the normalized final weight of the round up approach is 0.73, which is better than the 
final weight, 0.78, of the step 1 GA optimization. The result shows that the final design from step 1 was not the best 
design, and the solution had not fully converged yet. With step 2, the DOT optimization plus discretization 
approach, a better design was found with just an additional 297 function calls which is just about 8% of the total 
computational cost. The aeroelastic tailoring using the hybrid and discretization optimization approach studies is 
summarized in Table 7. The V-g and V-f plots of the final design are shown in Fig. 6. The final design has the same 
flutter mechanism, and the flutter speed constraint is satisfied. 

It was concluded that the hybrid and discretization optimization approach improves accuracy and computational 
efficiency of a global optimization algorithm. This study demonstrated that the effect of aeroelastic tailoring, which 
included ply orientation angles as design variables, produced an additional 10% weight reduction. A hybrid and 
discretization optimization approach improves accuracy and computational efficiency of a global optimization 
algorithm. 

B. Study 2: Flutter Mass Balancing 

As mentioned in Section III, the research objective of the X-56A aircraft is to demonstrate the simultaneous 
active suppression of three flutter mechanisms through the use of feedback controls. Thus, this aircraft is designed to 
have three flutter speeds within its flight envelope. The flutter speeds of the LMSW final design (non-validated 
model) with EFEW configuration were within the flight envelope as shown in Fig. 9. Unfortunately, after the ground 
vibration test (GVT) and subsequent model correlation (validated model), the second and third predicted flutter 
speeds of this validated model were too high or outside the flight envelope with the current propulsion systems and 
structural integrity. The normalized SWBT and AWBT flutter speeds at Mach 0.16 are about 1.48 and 1.68, 
respectively. The validated model was used as a baseline model in this flutter mass balancing study. Table 8 lists the 
baseline normalized flutter results at Mach 0.16 and the LMSW flutter requirements. The goal of this optimization 
study is to provide guidance to modify the wing design and move flutter speeds back into the flight envelope, so that 
the research objective can be accomplished. Since the vehicle is already fabricated, altering wing skin laminate (i.e. 
ply thickness and orientation angles) is not an option. Therefore, a mass balancing technique is employed. The 
ballast mass was modeled using concentrated mass elements in the finite element model. 

The optimization approach was chosen for the mass balancing study instead of trial and error. The objective 
function is to minimize the total ballast weight and/or the flutter speeds while satisfying flutter speed constraints. 
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Multi-points design is essential in this optimization since the total weight of the aircraft can be changed during 
flight. Two weight configurations, EFEW and FFFW, are taken into account in a single optimization run. The design 
requirements for the X-56A aircraft flutter speed and frequency constraints defined in Table 8 are used as inequality 
flutter constraints for the optimization. The design requirement for first flutter, SBFF, is 0.79 to 0.98, second flutter 
is 0.98 to 1.18, and third flutter is 0.98 to 1.30. Note that the predicted normalized flutter speeds of the baseline 
model configurations did not meet the flutter requirements. The V-g and V-f plots of the EFEW baseline model are 
shown in Fig. 8(a). The normalized SBFF flutter speed for EFEW is at approximately 1.13. The plots also show that 
coupling of the symmetric wing first bending (SW1B) and the symmetric wing first torsion (SW1T) mode has led to 
the SWBT normalized flutter speeds at approximately 1.48. The AWBT normalized flutter speeds occurred at 
approximately 1.68. 

To move flutter speed into the flight envelope using mass balancing optimization, DOT with CD Vs were used. 
Three design configurations and the ballast locations shown in Fig. 7 were investigated. To reduce the number of 
design variables and ensure symmetric design and mass distribution, design variable linking was used for left and 
right wings. The design variables selected in configuration 1 were six mass ballasts located primarily at the wing 
leading edge where each ranged up to 5 lb. For configuration 2, a total of thirteen design variables were selected, 
which included ten ballasts up to 5 lb each at the wing leading and trailing edges and one single ballast up to 20 lb at 
the nose of the center body. In configuration 3, the mass effectiveness was increased with a 25 -inch long aft wing tip 
boom and five lumped masses (ballast). The location of each ballast is defined within a five-inch segment, and each 
ballast can be moved with a one-inch increment. This trailing wing tip boom allows ballast to be added further 
behind the wing tip with a longer moment arm. The invention of the configuration 3 design was based on the 
optimization results of configurations 1 and 2. There were a total of eleven design variables in this configuration 
which included one center body nose ballast, five wing tip boom ballasts, and corresponding locations. The 
summary of each configuration is listed in Table 9. 

The final design variables after the optimization to reduce the flutter speed of configurations 1, 2, and 3 are listed 
in Tables 10, 11, and 12, respectively. The flutter results are listed in Table 13. Configuration 1 optimization results 
show that the design variable number 6 changed from 0.0 to 5.0, which represents a 5-lb ballast being added to the 
aft wing tip location. The third flutter speed was reduced from 1.68 to 1.59 and 1.67 for EFEW and FFFW cases, 
respectively, by the additional mass at aft wing tip. The optimization of configuration 1 and 2 has comparable results 
after adding 5 lb of ballast at the aft wing tip locations to reduce the third flutter speeds. These tendencies indicated 
that placing the mass toward the aft wing tip increases the mass effectiveness of the ballast. These results lead to the 
configuration 3 design, which is the trailing wing tip boom design configuration. The configuration 2 optimization 
results also show that adding 20-lb ballast at the nose of the center body can reduce the SBFF normalized speed 
from the configuration 1 values of 1.16 to 1.12 and 1.18 to 1.12 for EFEW and FFFW cases, respectively. 

The configuration 3 design allows for the placement of the mass ballast further behind the aft wing tip to 
increase the mass effectiveness. Since the goal is to move the second and third flutter speeds back into the flight 
region, the objective function is to minimize the first flutter speed. The constraints were the second and third flutter 
speeds defined in Table 8. The constraints were: 0.98 < V2 < 1.18, 0.98 < V3 < 1.3. All the flutter frequencies were 
used as inequality constraints as well. Three optimizations with different initial conditions were performed as shown 
in Table 12. In general, all three optimization cases achieve similar results when adding about 5-lb ballast to the 
most aft location of the trailing wing tip boom. Case 1 is the most practical design out of the three. The final designs 
of cases 1 and II are lighter than case III. Furthermore the case 1 final design is easier to implement than case II by 
adding just one ballast. Case I flutter results are listed in Table 13. The results show that the second and third flutter 
speeds can be moved into the desired speed region by using trailing wing tip boom ballast. Notice that for the EFEW 
configuration, the normalized SBFF speed is 1.13 and the normalized SWBT flutter speed is 1.11. The SWBT 
flutter speed was reduced too much by the 5-lb weight ballast to maintain the desired order of flutter mechanisms. 

As a result, a sensitivity study of the X-56A flutter characteristics using the trailing wing tip boom was 
performed. A 20-lb center body nose ballast was included in this study. The trailing wing tip boom tip ballast was 
varied from 1 to 5 lb with 1-lb weight increments. The flutter results are summarized in Table 14. The results 
indicated that to ensure the SBFF mechanism occurs before other flutter mechanisms, 3- or 4-lb instead of 5-lb 
ballast is recommended. With 5-lb ballast, the normalized EFEW AWBT flutter speed is 1.11 which is lower than 
the SBFF. The V-g and V-f plots of the EFEW model with 4-lb trailing wing tip boom ballast are shown in Fig. 
8(b). The second and third predicted flutter speeds are now within the flight envelope. Based on these results, either 
LMSW or NASA AFRC will design and implement the trailing wing tip boom with ballast weight on the flex wing 
in the near future for flight testing and perform active flutter suppression with all three flutter modes in the flight 
envelope. The non-dimensional flight envelope of the X-56A EFEW is shown in Fig. 9, which summarized the 
flutter results of the LMSW final design (non-validated), the baseline (validated), and flutter mass balancing. 
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VI. Conclusion 

An object-oriented MDAO tool that integrates aeroelastic effects has been developed and demonstrated. Two 
multidisciplinary design optimization studies were performed for the Multi Utility Technology Test-bed, X-56A 
aircraft. For all the sub-modules developed for this object-oriented MDAO tool, there is no restriction on the number 
of design variables and/or constraints as well as number of analytical configurations in a single optimization run. 
Although the sub-modules are developed for MSC Nastran and ZAERO, the sub-modules can be easily customized 
for other analytical tools. 

In this paper, the first example of the aeroelastic tailoring concept which exploited the anisotropic capabilities of 
fiber composite materials to minimize the structural weight while maintaining desired flutter speeds has been 
demonstrated using the X-56A aircraft final design for further optimization. Results indicated that the object- 
oriented MDAO tool with high fidelity structural analysis modules can be used for preliminary and detailed design 
tasks of aircraft. 

Results demonstrated that a genetic algorithm with discrete design variables is a beneficial approach for 
optimizing composite laminates with its ability to handle all types of design variables represented by realistic 
constraint sets in a finite element model, different material types, and manufacturing ability; thereby providing the 
flexibility needed to solve aeroelastic tailoring problems. Although a genetic algorithm is designed for global 
optimization, it is computationally inefficient and depends on the size of population and desired number of 
generations. In this study, a hybrid and discretization optimization approach is used to improve the computational 
efficiency problem with a global optimization. With the use of additional DOT optimization and discretization 
approaches following the genetic algorithm, the final design can be fine-tuned. 

The second example of the mass balancing optimization results showed that the X-56A symmetric body freedom 
flutter speed can be reduced by adding nose ballast. The second and third flutter speeds can be reduced to the desired 
region using a trailing wing tip boom. The results of this optimization study provided guidance to modify the wing 
design and alter the flutter speeds back into the flight envelope, so that the original research objective of the X-56A 
flight test (i.e. examine an active flutter suppression system to suppress all three flutter modes) can be accomplished 
successfully. This example also demonstrated that the object-oriented MDAO tool can handle multiple analytical 
configurations in a single optimization run. 
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Tables 


Table 1. Aeroelastic tailoring optimization constraints. 


Constraint 

Value 

Inequality Constraint 
Tolerance 

Comments 

Composite FI 1 

< 1.0 

0.001 

Based on 2.5g loads. Factor of safety of 1.5. 

Composite FI 2 

< 1.0 

0.001 

Based on -1.0 g loads. Factor of safety of 1.5. 

Buckling load factor 

> 1.0 

0.001 

Applied loads are based on 2.5 g loads. Factor 
of safety of 1.5. 

Normalized critical Flutter 
speed 

> 1.62* 

0.002 

At Mach 0.16 


*: Flutter speed constraint > 1.62: 

• Genetic algorithm: 162 V < 0.002; V > 1.617 

1.62 

• DOT: < (0.002 + 0.003) (CTMIN) ; V > 1.612 

1.62 


Table 2. Aeroelastic tailoring optimization cases (Constraints: V > 1.15V L ; BLF >1.0; and FI < 1.0). 


Case 

Design variable set 

Number of design 
variables 

Step 

Optimization 

descriptions 

Number of population and 
generation 

1 

Thickness 

12 

1 

GA + DDV 

200 x30 

2 

DOT + CDV 

N/A 

2 

Thickness + 
orientation 

24 

1 

GA + DDV 

200 x30 

2 

DOT + CDV 

N/A 


Table 3. Aeroelastic tailoring Case 1, ply thickness optimization case non-dimensional response values. 



Starting 

Step 1 

(GA +DDV) 

Step 2 

(DOT + CDV) 

Round down 

Round up 

Number of 
function calls 

N/A 

1406 

77 

N/A 

N/A 

Normalized 

weight 

reduction 

1.0 

0.88 

0.87 

0.86 

0.89 

Composite FI 1 

0.26 

0.28 

0.30 

0.30 

0.29 

Composite FI 2 

0.17 

0.19 

0.20 

0.20 

0.20 

Buckling load 
factor 

1.05 

1.14 

1.16 

1.18 

1.14 

Normalized 
critical flutter 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

1.69 

0.97 

1.62 

0.97 

1.61 

0.96 

1.60 

0.96 

1.64 

0.97 
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Table 4. Aeroelastic tailoring Case 2, ply thickness and orientation case non-dimensional response values. 



Starting 

Step 1 

(GA +DDV) 

Step 2 

(DOT + CDV) 

Round down 

Round up 

Number of 
function calls 

N/A 

3416 

297 

N/A 

N/A 

Normalized 

weight 

reduction 

1.0 

0.78 

0.72 

0.70 

0.73 

Composite FI 1 

0.26 

0.43 

0.40 

0.43 

0.39 

Composite FI 2 

0.087 

0.28 

0.27 

0.27 

0.26 

Buckling load 
factor 

1.05 

1.24 

1.26 

1.27 

1.25 

Normalized 
critical flutter 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

1.92 

0.83 

1.62 

0.57 

1.61 

0.56 

1.42 

0.64 

1.62 

0.63 


Table 5. Aeroelastic tailoring Case 1, ply thickness case design variable values (inch). 


Design 

variables 

Lower limit 

Starting 

After step 1 

After step 2 

Round 

down 

Round up 

Upper limit 

1 

0.0101 

0.1010 

0.0909 

0.0864 

0.0808 

0.0909 

0.1010 

2 

0.0101 

0.1010 

0.1010 

0.0945 

0.0909 

0.1010 

0.1010 

3 

0.0625 

0.6250 

0.1250 

0.0625 

0.0625 

0.0625 

0.6250 

4 

0.0101 

0.1010 

0.0707 

0.0537 

0.0505 

0.0606 

0.1010 

5 

0.0101 

0.1010 

0.0101 

0.0135 

0.0101 

0.0202 

0.1010 

6 

0.0101 

0.1010 

0.0808 

0.0630 

0.0606 

0.0707 

0.1010 

7 

0.0101 

0.1010 

0.0909 

0.1010 

0.1010 

0.1010 

0.1010 

8 

0.0101 

0.1010 

0.0808 

0.1010 

0.1010 

0.1010 

0.1010 

9 

0.0625 

0.6250 

0.6250 

0.6250 

0.6250 

0.6250 

0.6250 

10 

0.0101 

0.1010 

0.0606 

0.0510 

0.0505 

0.0606 

0.1010 

11 

0.0101 

0.1010 

0.0404 

0.0321 

0.0303 

0.0404 

0.1010 

12 

0.0101 

0.1010 

0.0707 

0.0606 

0.0606 

0.0606 

0.1010 
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Table 6. Aeroelastic tailoring Case 2, ply thickness and orientation case design variable values (thickness = 
inch and angle = degree). 


Design 

variables 

Lower limit 

Starting 

After step 1 

After step 2 

Round down 

Round up 

Upper limit 

1 

0.0101 

0.1010 

0.0606 

0.0349 

0.0303 

0.0404 

0.1010 

2 

0.0101 

0.1010 

0.0808 

0.0567 

0.0505 

0.0606 

0.1010 

3 

0.0625 

0.6250 

0.4375 

0.5177 

0.5000 

0.5625 

0.6250 

4 

-45.00 

45.00 

-45.00 

-45.00 

-45.00 

-45.00 

90.00 

5 

-45.00 

90.00 

-45.00 

-45.00 

-45.00 

-45.00 

90.00 

6 

-45.00 

45.00 

30.00 

31.84 

30.00 

30.00 

90.00 

7 

0.0101 

0.1010 

0.0505 

0.0167 

0.0101 

0.0202 

0.1010 

8 

0.0101 

0.1010 

0.0101 

0.0101 

0.0101 

0.0101 

0.1010 

9 

0.0101 

0.1010 

0.0101 

0.0101 

0.0101 

0.0101 

0.1010 

10 

-45.00 

45.00 

15.00 

36.31 

30.00 

30.00 

90.00 

11 

-45.00 

90.00 

0.00 

-19.39 

-15.00 

-15.00 

90.00 

12 

-45.00 

45.00 

90.00 

43.13 

45.00 

45.00 

90.00 

13 

0.0101 

0.1010 

0.0101 

0.0101 

0.0101 

0.0101 

0.1010 

14 

0.0101 

0.1010 

0.0808 

0.0803 

0.0707 

0.0808 

0.1010 

15 

0.0625 

0.6250 

0.5000 

0.4907 

0.4375 

0.5000 

0.6250 

16 

-45.00 

45.00 

90.00 

72.89 

75.00 

75.00 

90.00 

17 

-45.00 

90.00 

-45.00 

-45.00 

-45.00 

-45.00 

90.00 

18 

-45.00 

45.00 

75.00 

73.76 

75.00 

75.00 

90.00 

19 

0.0101 

0.1010 

0.0101 

0.0198 

0.0101 

0.0202 

0.1010 

20 

0.0101 

0.1010 

0.0303 

0.0181 

0.0101 

0.0202 

0.1010 

21 

0.0101 

0.1010 

0.0606 

0.0264 

0.0202 

0.0303 

0.1010 

22 

-45.00 

45.00 

-30.00 

-45.00 

-45.00 

-45.00 

90.00 

23 

-45.00 

90.00 

60.00 

44.03 

45.00 

45.00 

90.00 

24 

-45.00 

45.00 

-30.00 

-45.00 

-45.00 

-45.00 

90.00 


Table 7. Summary of the aeroelastic tailoring optimization study. 



Case 1, ply thickness 

Case 2, ply thickness and orientation 

Step 1 

(GA +DDV) 

After step 2 
round up 

Step 1 

(GA +DDV) 

After step 2 
round up 

Number of 
function calls 

1406 

77 

3416* 

297* 

Normalized 
weight reduction 

0.88** 

0.89 

0.78** 

0.73 

Composite FI 1 

0.28 

0.29 

0.43 

0.39 

Composite FI 2 

0.19 

0.20 

0.28 

0.26 

Buckling load 
factor 

1.14 

1.14 

1.24 

1.25 

Normalized 
critical flutter 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

Speed 

Freq. 

1.62 

0.97 

1.64 

0.97 

1.62 

0.57 

1.62 

0.63 


*: Effect of hybrid optimization + discretization: accelerate global optimizer; further improve design with about 


8% of additional function calls 

**: Effect of aeroelastic tailoring: 0.88 vs. 0.78 (10% more reduction by aeroelastic tailoring) 
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Table 8. Mass balancing baseline flutter predictions and requirements at Mach 0.16. 


Flutter 

mode 

Normalized speed 

Normalized frequency 

Lower 

bound 

requirement 

EFEW* 

FFFW* 

Upper 

bound 

requirement 

Lower 

bound 

requirement 

EFEW 

FFFW 

Upper 

bound 

requirement 

SBFF 

0.79 

1.13 

1.16 

0.98 

0.53 

0.68 

0.53 

1.76 

SWBT 

0.98 

1.48 

1.48 

1.18 

1.17 

2.34 

2.25 

2.35 

AWBT 

0.98 

1.68 

1.68 

1.30 

1.50 

1.52 

2.43 

3.52 


*: Flutter speed is beyond upper bound requirement 


Table 9. Flutter mass balancing design configurations descriptions. 


Configuration 

Number 
of design 
variables 

Objective 

Side constraints 

Comments 

Nose ballast (lb) 

Wing ballast (lb) 

Lower 

limit 

Upper 

limit 

Lower 

limit 

Upper 

limit 

1 

6 

Min. total 
ballast weight 
and target 
flutter speed 

N/A 

N/A 

0.0 

5.0 

Wing leading edge 
optimization: 5 wing leading 
edge and 1 aft wing tip 

2 

13 

Min. total 
ballast weight 
and target 
flutter speed 

0.0 

20.0 

0.0 

5.0 

Wing leading and trailing 
edge optimization: 1 nose, 6 
wing leading edge and 6 wing 
trailing edge locations 

3 

11 

Min. the first 
flutter speed. 

0.0 

20.0 

0.0 

5.0 

Aft wing tip boom 
optimization: 1 nose, 5 wing 
tip boom ballast and 5 wing 
tip ballast locations 


Table 10. Flutter mass balancing configuration 1 design variables. 


Configuration 1 

Design variable 

Initial value 

Final value 

Side constraints 

Lower limit 

Upper limit 

Wing leading edge ballast (lb) 

1 

0.0 

0.0 

0.0 

5.0 

2 

0.0 

0.0 

0.0 

5.0 

3 

0.0 

0.0 

0.0 

5.0 

4 

0.0 

0.0 

0.0 

5.0 

5 

0.0 

0.0 

0.0 

5.0 

Wing trailing edge ballast (lb) 

6 

0.0 

5.0 

0.0 

5.0 
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Table 11. Flutter mass balancing configuration 2 design variables. 



Configuration 2 

Design variable 

Initial value 

Final value 

Side constraints 

Lower limit 

Upper limit 


''Jose ballast (lb) 

1 

0.0 

20.0 

0.0 

20.0 

Wing leading edge ballast (lb) 

2 

0.0 

0.0 

0.0 

5.0 

3 

0.0 

0.0 

0.0 

5.0 

4 

0.0 

0.0 

0.0 

5.0 

5 

0.0 

0.0 

0.0 

5.0 

6 

0.0 

0.0 

0.0 

5.0 

7 

0.0 

0.0 

0.0 

5.0 

Wing trailing edge ballast (lb) 

8 

0.0 

0.0 

0.0 

5.0 

9 

0.0 

0.0 

0.0 

5.0 

10 

0.0 

0.0 

0.0 

5.0 

11 

0.0 

0.0 

0.0 

5.0 

12 

0.0 

0.0 

0.0 

5.0 

13 

0.0 

5.0 

0.0 

5.0 


Table 12. Flutter mass balancing configuration 3 design variables. 


Configuration 3 

Design 

variable 

Case I 

Case II 

Case III 

Side constraints 

Initial 

value 

Final 

value 

Initial 

value 

Final 

value 

Initial 

value 

Final 

value 

Lower 

limit 

Upper 

limit 

Nose ballast (lb 

) 

1 

20.0 

20.0 

20.0 

20.0 

20.0 

20.0 

0.0 

20.0 

Wing tip boom ballast (lb) 

2 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

5.0 

3 

0.0 

0.0 

4.8 

0.0 

1.0 

0.04 

0.0 

5.0 

4 

0.0 

0.0 

0.0 

0.0 

1.0 

0.04 

0.0 

5.0 

5 

0.0 

0.0 

0.0 

0.3 

1.0 

2.4 

0.0 

5.0 

6 

0.0 

5.0 

0.0 

4.7 

1.0 

5.0 

0.0 

5.0 

Wing tip ballast x location (inch) 

7 

216 

216 

216 

216 

216 

216 

211 

216 

8 

221 

221 

221 

221 

221 

221 

216 

221 

9 

226 

226 

226 

226 

226 

226 

221 

226 

10 

231 

231 

231 

231 

231 

231 

226 

231 

11 

236 

236 

236 

236 

236 

236 

231 

236 
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Table 13. Flutter mass balancing flutter results. 


Normalized flutter speeds 

Flutter mode 

Configuration 

1 

2 

3 - case I 

3 - case 11 

3 - case III 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

SBFF 

1.16 

1.18 

1.12 

1.12 

1.13 

1.14 

1.13 

1.14 

1.14 

1.14 

SWBT 

1.49 

1.67 

1.49 

1.67 

1.11* 

1.18 

1.12* 

1.18 

1.06* 

1.10* 

AWBT 

1.59 

1.57 

1.56 

1.55 

1.29 

1.26 

1.29 

1.26 

1.30 

1.26 

Normalized flutter frequency 

Flutter mode 

Configuration 

1 

2 

3 - case I 

3 - case 11 

3 - case III 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

EFEW 

FFFW 

SBFF 

0.66 

0.52 

0.71 

0.58 

0.72 

0.58 

0.72 

0.58 

0.71 

0.57 

SWBT 

1.28 

1.25 

1.28 

1.25 

1.07 

1.03 

1.07 

1.03 

0.99 

0.95 

AWBT 

2.09 

2.02 

2.07 

2.01 

1.57 

1.55 

1.57 

1.54 

1.44 

1.42 


*: Flutter happens before SBFF. 


Table 14. Wing tip boom flutter mass balancing sensitivity study flutter results (Note: includes 20 lb center 
body nose ballast). 


Wing tip boom ballast 
weight (lb) 

FFFW 

EFEW 

Flutter descriptions 

Normalized 

speed 

Normalized 

frequency 

Normalized 

speed 

Normalized 

frequency 

1 

1.12 

0.59 

1.11 

0.73 

SBFF 

1.37 

1.96 

1.38 

2.02 

AWBT 

1.60 

1.29 

1.39 

1.30 

SWBT 

2 

1.13 

0.59 

1.12 

0.73 

SBFF 

1.32 

1.82 

1.34 

1.87 

AWBT 

1.43 

1.19 

1.28 

1.23 

SWBT 

3 

1.13 

0.58 

1.12 

0.72 

SBFF 

1.29 

1.71 

1.31 

1.75 

AWBT 

1.32 

1.12 

1.21 

1.17 

SWBT 

4 

1.13 

0.58 

1.12 

0.72 

SBFF 

1.28 

1.62 

1.25 

1.66 

AWBT 

1.24 

1.07 

1.30 

1.12 

SWBT 

5 

1.14 

0.58 

1.13 

0.72 

SBFF 

1.26 

1.54 

1.11 

1.57 

AWBT 

1.18 

1.03 

1.29 

1.07 

SWBT 
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Figures 



Figure 1. Flowchart of the flutter analysis module in the object-oriented MDAO tool. 



Figure 2. Lockheed Martin Corporation X-56A aircraft (Figure courtesy: AFRL/Lockheed Martin). 
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Figure 3. X-56A baseline structural finite element model. 



Figure 4. X-56A baseline wing upper and lower skin design variable map. 
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Frequency, Hz Damping 



b) 




Non-dimensional speed 



0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Non-dimensionai speed 


Figure 5. The V-g and V-f plots for aeroelastic tailoring optimization initial design: (a) Case 1: ply thickness 
design variables case; (b) Case 2: ply thickness and orientation design variables cases. 
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Frequency, Hz Damping 


a) 


b) 






Figure 6. The V-g and V-f plots for aeroelastic tailoring optimization final design: (a) Case 1: ply thickness 
design variables case; (b) Case 2: ply thickness and orientation design variables case. 
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Frequency, Hz Damping 




c) Configuration 3, with a 25-in. railing wing tip boom (five 5-in. segments) 



Figure 7. Flutter mass balancing optimization design configurations. 


a) 


b) 






Non-dimensional speed 


Non-dimensional speed 


Figure 8. The V-g and V-f plots for empty fuel and empty water ballast configuration: (a) initial design; (b) 
final design with 20-lb nose and 4-lb trailing wing tip boom ballast. 
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Non-validated model 


• Symmetric body freedom 

flutter 

■ SWBT flutter 

♦ AWBT flutter 


Validated model 
O Symmetric body freedom 
flutter 

□ SWBT flutter 
O AWBT flutter 


After mass balancing 
O Symmetric body freedom 
flutter 

□ SWBT flutter 
O AWBT flutter 


0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Non-dimensional speed 


Figure 9. Non-dimensional X-56A empty fuel and empty water ballast configuration flight envelope. 
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